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Abstract. We define the Ladyzhenskaya-Lions exponent a L (n) = (2 + n)/4 
for Navier-Stokes equations with dissipation —(—A)" in 1", for all n > 2. We 
review the proof of strong global solvability when a > ct L (n), given smooth 
initial data. If the corresponding Euler equations for n > 2 were to allow 
uncontrolled growth of the enstrophy 5 II^ 7u IIl 2 ' * nen no globally controlled 
coercive quantity is currently known to exist that can regularize solutions of 
the Navier-Stokes equations for a < a L (n). The energy is critical under scale 
transformations only for a = a L (n). 



1. Introduction 

The Navier-Stokes equation for the velocity u of an incompressible fluid with 
viscous dissipation — (— A) a is given by 

fiii 

(1) — +u-Vu + \7p=-(~A) a u, 

at 

where u is a time-dependent divergence-free vector field in M™ with zero mean. 
The pressure p can be eliminated using Leray projections, since it merely serves to 
ensure the divergence free condition (see Equation (|20[) below). Given the initial 
condition 

(2) u(x,0)=u (x) 

where u (x) £ C ( ?°(R n ), the question is whether or not solutions remain smooth. 

Olga Ladyzhenskaya in the 1960s, proved global regularity for the two 
dimensional (n = 2) parabolic (a = 1) case. For n — 3 the challenge is to improve on 
the well known exponent a > 5/4 obtained long ago by Jacques-Louis Lions [IO]fTl]. 
still the best known result to date [H El [15]. One decade ago, Mattingly and 
Sinai |7J gave an elementary proof of this result and showed that there are no 
blow-ups and solutions remain smooth if a > 5/4. Later Katz and Pavlovic |6j 
extended the result by proving additional statements, concerning the Hausdorff 
dimension of the singular set at the time of first blow-up. Recently, in the context 
of theoretical physics, we have studied [T5] the problem in higher dimensions, in 
order to gain a better understanding of non-Newtonian turbulence and non-local 
anomalous diffusion of energy density in the Fourier domain in hyper-dissipative 
incompressible flows. 

We draw the reader's attention to how the lowest value of the exponent a known 
to ensure global regularity increases by 1/4 = 5/4— 1, as the number of dimensions 
increases from n = 2 to n = 3. The question that we address is whether or not 
this 1/4 increase per dimension is a more general principle with deeper significance. 
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Thus motivated, we define the Ladyzhenskaya-Lions exponent 




by linear extrapolation from n = 2 and n — 3. We then show how to generalize to 
all n the seminal regularity results of Ladyzhenskaya and Lions, starting with the 
following theorem: 

Theorem 1.1. If a > a L (n) then one has strong global solvability for the system 



In fact, we give two different proofs of this theorem. Although it may have gone 
somewhat unnoticed, Lions himself arrived at the result 40 years ago [11] (see also 
the paper by Guermond and Prudhomme [5]). Our own physical argument [T3] 
used L 1 estimates in the Fourier domain to bound the rate of growth of maximal 
norms of the velocity u and its spatial derivatives. Since then, we have been able to 
generalize to arbitrary n an elegant proof of Katz and Pavlovic [6] for the case n = 3, 
based on L 2 estimates. We present first the latter proof, because it uses standard 
methods of functional analysis. We then present the second proof, by reformulating 
in mathematical terms the physical arguments given in ref. |15] . This second proof 
is longer, but which nevertheless sheds further insight on how and why the cascades 
of energy (and vorticity) continue to bedevil attempts to prove or disprove global 
regularity for the important and well known special case a = 1 and n = 3. For 
higher dimensions n > 3 the nonlinearity can lead to even more violent cascade of 
energy, such that it becomes progressively harder to bound the uncontrolled growth 
of the enstrophy, defined as i|jVu||| 2 = — Au). 

The two methods of proof, though equivalent, shed light on different aspects of 
the problem and complement one another. Taken together, they seem to suggest 
that in higher dimensions n > 3 the Ladyzhenskaya-Lions exponent represents a 
genuine critical point [15] . Both methods rely ultimately (explicitly or implicitly) 
on conservation of energy to regularize the solutions. Solutions of the Navier-Stokes 
equations on a time interval [0,T] satisfy 



The second term on the right is called the dissipation term. The closely related 
Euler equations for inviscid flow do not have the dissipation term. 

We first remark on notation. We will use (, ) always to denote an L 2 inner 
product in space. We follow the notation used in ref. [6] throughout this paper for 
the expression A < B, to mean A < CB where C is some constant. This constant 
may depend on T, the norms of the initial values, and on an e which we keep fixed 
throughout this paper. 
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The following is a generalization to R™ of the known proof for n — 3. We briefly 
sketch the outline of the argument. By taking suitable inner products, one can 



m, 0) mR n . 




2. Proof based on L 2 



ESTIMATES 
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bound the rate of growth of LP Sobolev norms of u. Then, by invoking a carefully 
chosen Sobolev inequality, one can re-express these bounds in terms of L 2 norms, 
at which point the rest of the proof proceeds conventionally and one recovers global 
solvability. 

The dependence on n enters the final result solely through the Sobolev inequality. 
Most importantly, the dependence on n of a L (n) offers a new and different perspec- 
tive on the well known explanation of why the parabolic Navier-Stokes equations 
have global regularity in two but not in three space dimensions. Due to the close 
relationship between the Euler and Navier-Stokes equations, the literature does not 
always make clear the significant differences between them in two dimensions. One 
explanation argues that the for n = 2 vorticity is conserved, whereas for n = 3 vor- 
ticity and hence enstrophy can grow. This argument actually misses the point for 
the Navier-Stokes case, since in two dimensions the Euler equation itself has smooth 
solutions, rendering dissipation irrelevant. Plausible explanations have their basis 
on energy considerations. For n — 2 the direct energy cascade — if it counterfac- 
tually (i.e., hypothetically) existed in the Euler case — would remain weak enough 
that dissipation with a = 1 can control it. But for n > 2 the energy can cascade 
more violently, such that dissipation with a = 1 cannot necessarily control it. 

In this context, the following proof provides a clear interpretation, viz. that 
Sobolev embedding becomes more "costly" for larger n. By costly, we mean that 
the number of derivatives necessary for the embedding increases with n. Here, 
energy enters the picture (via the norm ||u|j^ 2 ), and the role of n becomes clear. 

Proof of Theorem ] Let H° and W"' p denote the L 2 and LP Sobolev spaces, 
respectively, over R" with [3 derivatives. Taking the inner product by pairing 
Equation (jTJ) with u, we get 

ld_ 

Hence, if the solution u remains smooth up to time T, we have the estimate 
(5) [ T \\uf Ha dt<{l + T). 



(4) 7^IMI| 2 = -II(-A)* u\\ 2 L2 < -\\u\\ 2 Ha + \\u\\ 2 L2 



d( 1 1 u II ) 

We next pair Equation ([T]) with (— A)^u in order to estimate dt Hf:> . We obtain 

(6) ^H(- A ) lu ll" 2 + < M ' Vu ' (- A )" u > = 
We next estimate the nonlinear term 

(7) (u ■ Vu, (-A)^) = ((-A)4(« • Vu), (-A)^u). 

Recall that u is divergence free. We can bound the absolute value of the left hand 
side of Equation (0 by 

(8) ||u • \7u\\ H n\\u\\ H u , 

using the Cauchy-Schwarz inequality. Next, we can apply the following well known 
generalization 1 of Holder's inequality: 



I/1/2 Hi* < II/iIUm ll/a 



1 _ 1 1 

P Pi Pi 



So we can now bound © by 

(9) ||u|| w ra,p||u|| 1 vfl+i, 9 ||u|| Hi B 
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or equivalently, following ref [6] , by 

(10) IMIw i .*IHIw''.«IMI.ff0 i 

with 

111 

2 p q 

Now comes the crucial step, where n enters the picture. We twice invoke the 
Sobolev embedding theorem to obtain, 

(12) \\u\\ w a, P \\u\\ w e+i,4u\\ HIJ <\\u\\ H «\\u\\ Ha +p\\u\\ Hf i , 

where, using (fTTj) . the embedding is seen to be completely continuous provided that 

(13) 2a- 1 > n/2 . 

The embedding "costs" n/2 derivatives whereas we can "spend" 2a— 1 derivatives. 
Remarkably, this is the only step of the argument involving n. We note that the 
same bound is obtained by using (fT0|) instead of J9} in (fi"2|) . 

The rest of the argument proceeds straightforwardly. From Cauchy-Schwarz, we 
get, 

\\u\\ H 4u\\ Ha+ 4u\\ H , < S\\u\\ 2 Ha+ , + j\\u\\ 2 Ha \\u\\ 2 Hft . 
Combining this with Equation we get 

j t \\u\\^<\\u\\U\u\\%, + \\u\\h. 

Together with © and invoking Gronwall's inequality gives global solvability. The 
sole sufficient condition ([TBI can be rewritten as, 



(14) a > — 



□ 



3. Proof based on L 1 (or L°°) estimates 

We next give a second proof corresponding to the original argument presented 
in our earlier study, ref. [TS]. The latter argument has its basis on the analogy 
or connection [15] between the Navier-Stokes equation and nonlocal Fokker-Planck 
equations for describing anomalous diffusion of ensembles of random walkers, such 
as found in Levy flights [12]. Fourier transforming the Euler equation converts 
the nonlinear term into a non-local one. Taking an inner product of the Fourier 
transformed equation with the Fourier transform of u, one obtains an equation 
which conserves energy, in a manner analogous to how Fokker-Planck equations 
conserve probability. A well known fact about fractional or nonlocal Fokker-Planck 
equations from the study of Levy flights concerns the relationship between fat tails 
in the probability density function and smoothness of the characteristic functions: 
the characteristic function may not be smooth for fat tailed probability density 
functions having diverging moments. Dissipation in the Navier-Stokes equations, 
of course, makes kinetic energy decay and adds an element of exponential damping 
to the Fourier transformed equation. We can use this damping to control and 
bound the growth of suitable norms. For convenience we chose to study the L 1 
norms involving the Fourier transform of u. The bounds then allow us to obtain 
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L°° estimates for u and its derivatives. Hence, we actually work in the classical 
Holder space rather than in Sobolev space. 

We first briefly sketch the outline of the proof. We take the Fourier transform of 
the Navier-Stokes equation and then bound the rate of growth of the Fourier trans- 
form u(k 7 .) of u(x, .), where k is the Fourier conjugate of x. Using this inequality, 
we then bound the "moments" of the absolute value of u. These "moments" are 
defined in a manner analogously to statistical moments of probability density func- 
tions. Conservation of energy plays an explicit and important role in the argument, 
because it leads to estimates of lower order derivatives of u in terms of higher order 
derivatives, which leads to the final result. 

One aspect of the problem which becomes more clear in this argument concerns 
the crucial role played by scale invariance. We will show in the concluding remarks 
that the critical value a L (n) = (2 + n)/4 for the marginal case corresponds to the 
scenario where the dissipation effects and the nonlinear effects (or nonlocal effects 
in the Fourier domain) have the same scaling. The number of dimensions n enters 
the picture in the step involving conservation of energy. This longer proof provides 
additional and complementary insight into how and why conservation of energy 
has different effects for different n. (In fact, both the divergence-free condition 
of incompressibility and conservation of energy become less relevant as n tends to 
infinity. ) 

In what follows Ui denotes the Cartesian component i — 1, 2, . . . n of u and \k\ is 
the length of the n-dimensional vector k. 

Proposition 3.1. The Fourier transform Ui(k, .) of Ui(x, .) satisfies 

llli^lL < EEfTjIli^ll^lli^r-^ll^ 

3 t 

(15) -\\\k\ 2a+m ,u t \\ L1 +C hm (t) . 

for any integer m > 0, where Ci tTn depend only on pressure p. 

Proof of Proposition \S.1[ Fourier transforming Equation ([1]) to eliminate the space 
derivatives, we obtain 

d d 

(16) — Ui(k,t) +E"i(M) * (ify«i(M)) = -|£f Q u*(M) ~ ikip(k,t) 

3 

with the nonlinear term becoming a nonlocal convolution term in Fourier space. 
Here i = \f— 1, whereas i denotes a component index. We now note that for any 
complex function g (i), if dg/dt = A(t) - B{t)g(t) and B > then f t \g\ < \A\ - B\g\. 
Identifying g with Uj and B with |/c| 2a , we get from Equation (|I6I) . 

- |fc| 2a |fii(M)| + |**p(M)| 



(17) -\ui(k,t)\ < 



Uj (fc, t) * (ikjUi(k, i)) 
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Multiplying by \k\ m (to = 0, 1, 2, . . .) and integrating out k over the entire Fourier 
domain R™ , we get, 



^ f dk ifcn^(M)i < f dk ifcrv 



(18) 
Hence, 



+ 



v,j(k,t) * (kjUi(k,t)j 

dk \k\ 2a+m \ui{k,t)\ 
dk |fc| m (|fcip(M)|) • 



— I 

dV 



*»ll£i 



< E/ dk '\W k + k '\ m ^\\^ \k'jUi{k',t)\ 



\\k\ 2a + m Vi\\ L , + C itm (t) 



/m 
d^ 
i) 



where Ci <m (t) 

d 
di 

(19) 



\\\k\%\\ L1 \k'\ m -*\k'M k ',t) 
-|||*| 201+m 5<|| L1 + C iini (t) 
|fc| m+1 p|| T 1 . The final integration leads to 



[i*r*L ^ EE 



I II rill 1*1 



Mlii 



\\k\ 2a+m ,u4 Ll +C hm (t) 



□ 



The pressure satisfies an elliptic Poisson equation, which we can invert to obtain, 

(20) p = -A^Tr (Vu) 2 . 

It corresponds to the purely non-divergence-free component of the nonlinearity. 
Henceforth, we ignore Ci >m since pressure plays no relevant role, as discussed earlier 
in the context of Leray projections. 

Proposition 13.11 allows us to control the growth of certain norms. To prove 
the existence of smooth solutions, we need the term with 2a in the exponent to 
dominate over the sums of product term (i.e., the inertial or nonlinear term). If 
this condition is satisfied, then the norms remain bounded. We next develop a few 
useful lemmas, which we prove for completeness. 

Lemma 3.2. Let f denote the Fourier transform of a function f : R" — ► R n . If 

e + n/2 

feC m ,then \k\ e f < \k\ m f mU forl<m. . 

L 1 L 1 

Proof. Scale transformations by a factor A show that the inner product 

(21) </,/) = (/»,/») 

(22) = (A-«/ 2 />/A),A-"/ 2 />/A)) , 
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so that f(k) and f\(k) = A n / 2 f(k/\) have the same inner product in L 2 = H°. 
Now 



(23) 
(24) 



\Kfx 
\k\ m fx 



L 1 



So for sufficiently large A we get 
1*1'/ 



\k\ £ h 



< 



L 1 



|fc|"7 

|fc| m /A 



L 1 



Moreover 



(25) 



l/(t+n/2) 
L 1 



\Wfx 



L 1 

l/(i+n/2) 



\k\ m f 



l/(m+n/2) 



from which it follows that || |/c|^/ < C 
the lemma. 



i*r/ 



l/(m+n/2) 



L 1 



for some finite C, proving 



□ 



Lemma 3.3. Let [3 = (f3\, /3%, ■ ■ ■ , Pn) be a multi-index and let \/3\ denote the sum of 



the components, according to convention. If 
Proof. 

( 27r )«/2|/(/5)( x ) I < 



\kf\f{k) 



< +oo then f £ C^l 



L 1 



(26) 



(2^)"/ 2 ||/^|| 



< 



|/»| le^ldfc 

\kf r f(k) 



L 1 



Hence, / belongs to the Sobolev space M^' /3 ''°°, consequently to the classical Holder 
space 

□ 



Proof of Theorem ll.il We first use Lemma 13.21 to bound the norms in the inertial 
terms: 

l+n/2 

(27) \\\k\%\\ T , < \k\ 



2a + m + n/2 

L 1 



(28) 



im-f+l- 



Uj 



< 



\k\ 



-i+l + n/2 



2a + m + n/2 



L 1 



From Proposition 13. li we retain control so long as the sum of the exponents on 
the right hand sides of (|27|) and (|28[) is less than unity. The condition to bound the 
rate of growth of the relevant L 1 norms then becomes 

(29) W + 1 + " < 1 

v ' 2a + m + n/2 

from which we recover the condition (|14[) . Lemma |3~31 guarantees global solvability, 
completing the proof. 

□ 



8 



T. M. VISWANATHAN AND G. M. VISWANATHAN 



4. Concluding remarks 



Attention is traditionally given in the case n = 3 to vorticity. In our proofs the 
key quantity of interest in the Navier-Stokes system clearly is the energy. As noted 
by Tao [T3l [14], the energy is a significant globally controlled coercive quantity 
which for a = 1 is critical for n = 2 but supercritical for n > 2. However, we note 
that for a > a h (n), the energy remains both coercive and either critical (a = a L (n)) 
or subcritical (a > a L (n)). We can state these ideas formally: 

Theorem 4.1. The kinetic energy ^\\u(.,T)\\ 2 L2 is critical under scale transforma- 
tions only for a — a L (n). 

Proof. If u is a solution, then 



We thus see that the value a L (n) = (2 + n)/4 corresponds to a critical value, for 
which the energy remains invariant under scale transformations, i.e. E\ = E. □ 

Remark 4.2. We briefly comment on the marginal case a = a L . We know from 
Proposition 13.11 and Equations (|2T[) and ([28]) that the problem for the marginal 
case arises because the nonlinear and dissipation terms scale identically and thus 
balance each other, so that it is not immediately clear which term dominates. Any 
small additional consideration can tip the balance. To address this issue, we recall 
that 



for T2 > T\ 1 i.e. the kinetic energy is a strictly decreasing function of time. This 
reduction in kinetic energy, insignificant for the non-marginal cases a < a h (n), can 
by itself rule out any blow-up for the marginal case, as suggested by the criticality 
condition in Equation (|31[) . 

We could, alternatively, describe the situation from a different perspective, viz., 
blow-ups are energetically forbidden for a > a L (n), but not for a < a L (n). The 
question is what happens exactly at the critical point. In principle, if the kinetic 
energy were a constant of the motion, then we would not be able to rule out the 
possibility that almost all the energy might become concentrated into ever smaller 
regions of space and at an increasingly faster rate. In this scenario, we would not 
be able to rule out finite time singularities. However, kinetic energy is a strictly 
decreasing function of time, hence in the time that it takes to concentrate the energy 
into a smaller region, enough energy will have dissipated to prevent the formation 
of singularities in finite time, because not only do dissipation and nonlinearity 
have the same scaling for a = a L (n), but in addition the kinetic energy decays 
asymptotically to zero at large times. 

Finally, we comment on the hypothetical worst-case scenario where the corre- 
sponding Euler equations for inviscid flows lead to blow-ups. No globally controlled 
coercive critical or subcritical quantity is known to exist for a < a L , therefore we 
are forced to consider the possibility that some initial conditions may very well lead 
to blow-ups in the Navier-Stokes case for a < a L in this hypothetical scenario. 



(30) 




is also a solution [15] . The energy E\ of the new solution scales to 
(31) E x = \ 2 ~72^T)E. 



(32) 



IK,T a )||i a <||«(.,T 1 )||i 
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For n — 3, the Euler equations may allow for the uncontrolled growth of the 
vorticity u> — V x u due to the vorticity stretching mechanism [4j [2j EJ - The 
enstrophy 

(33) i|M|| 2 = i( V x u,u>) = ^<u,V xu) = --(u,Au) = -\\Vu\\ L , 
can thus conceivably blow up because it also is not conserved: 

(34) ^IMl£* = (w-Vu,a;) ■ 

If it could be shown that a blow-up is possible for the Euler equation, then there 
may indeed be no way to regularize this worst-case scenario with a dissipation 
having a < a L (3). In this context, we recall the statement of Peter Constantin [3]: 
"It is no exaggeration to say that the Euler equations are the very core of fluid 
dynamics." 

The same argument can be generalized to n > 3: if |jVw|j^ 2 blows up in the 
Euler case, no way is (currently) known to regularize solutions via dissipation for 
a < a L (n). Although the question remains open, we speculate that the || Vm||| 2 can 
blow up in solutions of the Euler equations in three and higher dimensions. This 
intuition leads to our final statement: 

Conjecture 4.3. The Lady zhenskay a- Lions exponent a L (n) is critical in the sense 
that it separates two regions. Solutions of the hyper- dissipative Navier-Stokes equa- 
tions in dimensions n > 3 remain smooth for a > cc L (n), whereas for any a < a L (n), 
finite time singularities of ||Vu||| 2 are possible. 
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